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We investigate accreting disk systems with polytropic gas in Keplerian 
motion. Numerical data and partial analytic results show that the self- 
gravitation of the disk speeds up its rotation - its rotational frequency is 
larger than that given by the well known strictly Keplerian formula that 
takes into account the central mass only. Thus determination of central 
mass in systems with massive disks requires great care - the strictly Keple- 
rian formula yields only an upper bound. The effect of self-gravity depends 
on geometric aspects of disk configurations. Disk systems with a small (circa 
10 -4 ) ratio of the innermost radius to the outermost disk radius have the 
central mass close to the upper limit, but if this ratio is of the order of 
unity then the central mass can be smaller by many orders of magnitude 
from this bound. 

1. Introduction 

The literature on those accretion disk systems where the gas exhibits 
Keplerian motion ui 2 = uj^/r 3 (here to is the angular velocity of the gas, 
r denotes the cylindrical radius, and ujq is a constant) brings a confusing 
message. Usually the velocity of rotation of such disks is understood to be 
dominated by a heavy central mass M c , that is the strictly Keplerian equal- 
ity M c = Uq/G, where G is the gravitational constant, holds true. The disk 
is assumed to be light. The analysis of the massive black hole in the mas- 
ing nucleus of NGC 4258 exemplifies this interpretation [1] (Interestingly, 
modelling of this system confirms that the disk is light [2]). On the other 
hand, there are sporadic observations in the literature that even systems 
with massive disks can conform to the Keplerian rotation curve [10]. In 
such a case the central mass cannot be deduced solely from the Keplerian 
frequency. Quite recently Hure et al. [11] investigated thin dust disks and, 
assuming a monomial density profile, arrived to a conclusion that the rota- 
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tional Keplerian velocity gives only information on the enclosed mass of the 
system, i.e., a sum of the central mass M c and the mass of the disk 

In this paper we deal with self-gravitating barotropic disks in Keplerian 
motion. We show that, contrary to the folk belief, the Keplerian rotation 
curve in pure hydrodynamics does not exclude heavy disks. The angular 
frequency of the rotation curve does depend on the two masses, but they do 
not contribute additively. We start the discussion with analytic results con- 
cerning configurations with strictly (or almost strictly) Keplerian rotation 
law, that is where loq is equal (or close) to GM C . This happens for elongated 
disks, for which the ratio of the innermost radius r; n to the outermost one 
r ou t is small (or tends to zero), or in the test-fluid regime, that is when the 
gravity of the disk can be neglected. 

In numerical calculations we specialise to polytropic equations of state. 
The results show that the angular frequency wo satisfies GM C < ujq < 
G(M C + M(j), but the precise value of uq depends on the geometry of the 
disk, in particular on the ratio ri n /r out . If r- m /r ou t is small enough, then the 
central mass is reasonably well approximated by the Keplerian expression 
U)q/G. A significant deviation from this value can occur for elongated disks 
when their mass exceeds the central mass by a few orders of magnitude. In 
contrast to that, when r- m /r out = 1/2, the ratio ujq/G yields a value close 
to the central mass M c only when in addition <C M c ; in such 
the disk must be geometrically thin. In this specific example, obtained for 
the polytropic index T = 5/3, if ~ M c , the disk becomes thick and 
ojq/(GM c ) ~ 1.3 instead of 1; the latter ratio grows rapidly with the further 
increase of 

Models of equilibrium toroidal or disk-like figures constitute a great sim- 
plification of the astrophysical reality. While such configurations appear in 
the late stages of current simulations of binary mergers, the theory of thick 
accretion disks would require more complex modelling. Nevertheless, we are 
convinced that a good understanding of the properties of self-gravitating 
equilibrium disk-like figures is a prerequisite for the construction of more 
realistic models (cf. [7]). We return to this and connected issues in Sec. [6} 

2. Equations 

We consider a stationary, self-gravitating disk of gas rotating around a 
central point mass M c . The gravitational potential of the system can be 
written as the superposition <3? = — GM c /\x\ + $ g , where the potential <&g 
is due to the gravity of the disk. It satisfies the Poisson equation 



A$ g = AirGp, 
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where p denotes the mass density of the gas. Here x = (x, y, z) are the Carte- 
sian coordinates originating at the central point-like mass. In the following 
we will also work in the cylindrical coordinates (r, <J), z). 

Assume that the disk is axially symmetric and occupies a finite volume. 
The last demand imposes a condition onto equations of state, as discussed in 
Appendix A. For technical reasons we will also assume that the fluid region 
is connected, and that it is equatorially symmetric. The velocity of the gas 
is U = u(r, z)d ( f ) or equivalently U = roj(r, z)e^, where denotes the unit 
azimuthal vector. A classical result due to Poincare and Wavre states that 
for a barotropic equation of state cj depends only on the cylindrical radius 
r [T7]. We will deal later with polytropic equations of state p = Kp r , where 
p is the gas pressure, and K and V are constant. The Euler equations 

Vp + p(U • V)U + pV$ = 

can be integrated, yielding 

h + $ c + $ = C (2) 

in the closure of a region where p is nonzero . Here h denotes the specific 
enthalpy of the fluid: dh = dp/ p, and 

/r 
dr'r'oj 2 (r') 

is the centrifugal potential. The structure of the disk can be obtained from 
Eqs. ([!]) and Q provided that the equation of state and the rotation law 
u = oj{r) are known. We would like to point out that the process of solving 
ofEq. g consists in the simultaneous finding of the unknown functions and 
of the volume where the enthalpy h is nonnegative. 



3. Keplerian rotation laws 

In the rest of this paper we investigate systems with the Keplerian ro- 
tation laws 



They are of special interest, because they agree with many observations. 
Almost all known stellar black holes with gaseous disks and a few of AGN's 
obey this law. 

Suppose now that in an accretion system the motion of a gaseous disk 
conforms to a Keplerian rotation law. The mass of gas (usually unknown 
to an observer) can impact the rotation curve. How legitimate is the ( often 
made) guess that the motion is almost strictly Keplerian, ujq ~ GM C , that 



4 



MMP printed on November 26, 2012 



allows one to infer the central mass from the observed rotation law? Analytic 
answers can be given in a few special cases, i) A thin gaseous Keplerian disk, 
with a mass negligibly small in comparison to the central mass M c , exhibits 
strictly Keplerian motion, i.e., loq = GM C . This is a well known fact, and 
we give in Sec. 3.1 a proof, only for the sake of completeness. Thus, in this 
case, one can obtain the value of the central mass simply by measuring 
the rotation frequency, ii) For thick and heavy disks the situation is more 
complex. There is still a possibility to prove analytically the approximate 
equality ojq ~ GM C , but only in special cases and under certain assumptions 
on disk's geometry. These results are discussed in Sec. 3.2 and 3.3. 

In the generic case the presence of massive disks influences the rotation 
curve, or strictly saying, the frequency parameter. Our numerical analysis, 
reported in later sections, suggests that the difference luq — GM C is strictly 
positive and increases with the increase of Md/M c , where = f d?xp is 
the mass of the disk. A result of Hure et al. [Tl] (obtained for dust) suggests 
that the difference ujq — GM C should be proportional to the ratio M^/M c . 
We show later that this is not true for polytropes. 

The centrifugal potential for the rotation law ^ can be written in the 
form $ c = iOy/r. Suppose that the disk is finite and extends from the inner- 
most cylinder labelled by r; n to the outermost cylinder at r = r out , and the 
innermost and outermost points of the disk are located on the equatorial 
plane. It follows from Eq. ^ that 

u$ = GM C - r in $ g (r in ) + r- m C = GM C - r out $g(r out ) + r out C, (4) 

or 

.o^GM^ ^^-^^l (5) 

Tout ?*in 

because h should vanish both for r- m and r ou t. Here <3? g (ri n ) = <3? g (r = n n , z = 
0) and $ g (r out ) = $ g (r = r out , z = 0). 

There are simple analytic arguments - we discuss them in Sec. 3.4 - 
showing that Wg exceeds GM C: when the constant C > 0. In the most of 
this paper, however, we deal with finite-volume configurations for which C 
can be negative, and the analytic argument of Sec. 3.4 does not work (the 
relation between the equation of state, Keplerian rotation, and the sign of C 
is discussed in Appendix A). Nevertheless, our numerical data reported in 
later sections do suggest that in fact GM C < ujq < G(M C +Md). We state this 
as a conjecture: Barotropic Keplerian disks rotate with the angular velocity 
exceeding (GMc/r 3 ) 1 / 2 , but smaller than (G(M C + M d ) /r 3 ) 1 / 2 . 

3.1. Test- fluid limit 

Massless Keplerian disks must be i) strictly Keplerian and ii) infinitely 
thin. 
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Indeed, choose disk configurations with fixed r; n and r out , but let the 
density p ma , x — > 0. In this limit, the potential $ g also vanishes (<6 g satisfies 
Eq. (I), and it is normalised to zero at infinity). Part i) follows now from 
Eq. (5), which yields u)q — > GM C . 

To prove part ii), let us observe that Eq. Q yields C = in the limit 
of $ g — > 0. In this case Eq. (pi) should be written as 



h + GM C f- — J—^ ) =0, 

and this form admits (notice that h — > when p max 0) only infinitely 
thin disks at z = 0. 



3.2. Disks with small inner radii 

Let us fix the maximum density within the disk, say p = p max , and let 
r out be also fixed. This choice is justified by the scaling symmetry of the 
equations (see Sec. [4]). Then, it follows from Eq. ^ that ujq — > GM C as 
Hn - > 0. To show this, it suffices to note that is a bounded function on 
M . To be more precise: |^ , g(T ou t) — < ^g(nn)| must be finite when /3 max is 
fixed. 

The above result can be intuitively understood, because for r small 
enough, the gravitational potential is always dominated by the divergent 
term — GM c /|x|. In practice, if the ratio ri n /r ou t is sufficiently small, then 
the rotation curve of the gaseous disk should be influenced "almost exclu- 
sively" by the central mass. We show later, solving specific examples, that 
this is true. 



3.3. Ring-like disks 

One can also show that ojq — > GM C as r- m — > r out , that is when disk 
becomes ring-like with a small inner radius (r ou t— r- m ) /2. The proof proceeds 
according to the squeeze theorem. For the sake of simplicity the second part 
of the reasoning (the bound from below) will be presented for polytropic 
equations of state, but the proof can be reformulated for general barotropes 
satisfying certain additional conditions. 

The starting point is the virial relation formulated in Mach |13j for the 
system consisting of a steady disk and a point mass. It yields 

i J p$ g d 3 x - J p ( ^^d 3 x + J p\U\ 2 d 3 x + 3 J pd 3 x = 0. (6) 

For Keplerian rotation this equation can be written as 

\ fp%d 3 x+ L(^-9^A d 3 x + 3 fpd 3 x = 0. (7) 



G 
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Now, because p > we have 
and, since |x| > r, 

(wg - GM C ) y ^ 3 x < - X - J p$ g d 3 x. 

Let $ gi min denote the minimum value of the potential <& g in the region where 
p ^ 0. We have 



w n 2 - GM C 



Tout 



J pd 3 x < ~® gtmia J pd 3 x, 



and thus 



2 1 

UJ - GM C < --$g, m i n rout- 



For r; n —7- r ou t the mass of the disk tends to zero (note that we keep /9 max 
fixed, and the volume of the disk tends to zero), and so does the gravitational 
potential <J> g , m ax- Thus, in the limit we have ujq — GM C < 0. 

On the other hand, the same formulation of the virial theorem gives 

J J \ T |x| J \T'va ?"out / J 

which, for polytropic equations of state implies that 

-3^<"t_^. 

Hn Tout 

From the integrated Euler equation we have 

KT r l * ^ w o GM c _ ^ , , w 2 GM C 
+ *.M + f - — = C = * 5 (r ou ,) + J- - — , 

where 3> g (r) = 3> g (r = r, z = 0), and f denotes a cylindrical radius of 
the point where the density reaches its maximum. For simplicity we have 
assumed that the maximum lies on the equatorial plane. The second equality 
expresses the Euler equation at the outer boundary. Now, for r m — > r out 
we have f — > r out as well. In this limit one obtains K = 0, so that < 
Wq — GM C < 0. This concludes the proof. 
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3-4- Infinite disks 

If disks are infinitely extended and posses a finite mass, then C = 0. (An 
infinite extension imposes severe restrictions onto the allowed equations of 
state - see also the discussion in Appendix A.) We claim that in this case 
ool > GM C . 

Indeed, Eq. ([I]) yields, from the maximum principle [8], that the potential 
<£ g is nonpositive everywhere and vanishes at spatial infinity. Thus Eq. Q 
implies that the rotational angular frequency must exceed the Keplerian 
value and loq > GM C . 

The same argument would work for the constant C > 0, but it is not 
clear whether such configurations exist. 

4. Rescaling of equations 

It is convenient to transform Eqs. ([I]) and ^ into dimensionless forms. 
The quantity u = Gr^ ut p ma , x has the dimension of potentials. We define 
rescaled dimensionless potentials i& c = ^ c /u, $ = $/"U, <£ g = 5> g /«. The 
scaled enthalpy is given by h = h/u while the transformed density reads 
p = p/pmax- New spatial coordinates are defined as x = x/r out . Introducing 
M c = GM C / (pmax^out)) we can s P n * the new potentials as $ = — M c /|x|+$ g . 
A similar trick has been done in [3] and [9]. The disk-related gravitational 
potential $ g satisfies 

M g = 4^5, (8) 

where A is the laplacian with respect to the new coordinates x. The Euler 
equation Q reads 

h + <l c + $ = C. (9) 

For the polytropic equations of state, the enthalpy h can be expressed as 
h = (KT/(T - l))p r_1 , and h = (^r/(T - l))p r -\ where constants K 
and K are related by K = Kp^^/u. Finally, for the Keplerian rotation 
law we can define $ c = Cj^/f. In this rescaled notation the Keplerian case 
loq = GM C corresponds to cDq = M c . 

In the new variables, the set of parameters specifying the solution is 
reduced to T, M c and f- m = r- u Jr OU f Taking large values of M c corresponds 
to test-fluid solutions. Small values of M c yield solutions with the massive 
disk. 

5. Numerical examples 

The numerical method of this paper follows the classic Self-Consistent 
Field (SCF) scheme jTU 0]. Eqs. ^ and ([£]) are solved iteratively on 
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a two dimensional grid. In every iteration step the gravitational potential 
of the disk is found from Eq. ([8]), assuming that the density is known. A 
new density distribution is then obtained directly from Eq. ([9]). Eq. Q is 
solved by the Green function formula, where the integral kernel is expanded 
in Legendre polynomials. In this approach the boundary conditions for the 
potential are automatically satisfied, and one can restrict the numerical 
grid to a region containing the disk. We would like to point out that this 
method requires renormalisation of the constant C appearing in @ and the 
polytropic constant K. They are adjusted in each iteration step in order 
to agree with the assumed parameters of the solution - r- m , r ou t and the 
maximal value of p equal to 1. The constant u)o is computed from ^ at 
each iteration step. 

There exist modern numerical methods that are characterised by bet- 
ter convergence properties [61 13] . Their implementation would probably be 
necessary for non-Keplerian rotation laws. For Keplerian rotation laws nu- 
merical solutions converge well even with the SCF method - one can obtain 
a solution for nearly any set of parameters. The drawback of the method is 
the need to employ a large number of Legendre polynomials for thin disks. 
One is able, using a parallel version of the code and optimisations described 
in [12] , to find accurate solutions using up to 400 of Legendre polynomials 
and large grids (up to 5000 x 5000). Such resolutions are permitted because 
of the simplicity of the SCF method. 

Our numerical solutions successfully passed a virial test in the new for- 
mulation that includes a point mass |13j . The virial theorem ^ can be 
written as 

■Spot + 2-Ekin + 2£'therm = 0, 

where E pot = J d 3 xp$ s /2-GM c J d 3 xp/\x\, E kin = J d 3 xp\XJ\ 2 /2, E therm = 
J d 3 x3p/2. For the tested quantity we choose 

£v = | -Spot + 2-Ekin + 2i?therm I / 1 -^pot | • 

We have obtained a satisfactory accuracy, documented by the smallness of 
the ratio e v . The value of e v depends on the resolution of the numerical grid. 
Appendix B shows that e v decreases from the value of 10 -4 to 10 -8 with 
the finessing of the resolution. 

We performed a few thousands runs for polytropic disks, with two dif- 
ferent polytropic indices T = 5/3 and T = 4/3. The results that we would 
like to discuss first are presented on Figs. [I}j6| Figures [T] and [2] show the 
graphs of Uq — M c and Md versus rj n for a couple of different values of M c . 
They confirm the conjecture stated in Sec. [3j which in the rescaled variables 
can be written as M c < u% < M c + Md or < (in — M c < Md. It is also 



clear that ujq — > M c for f- m — > and f; n — > 1, as it was proven in Sec. 3.2 
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Fig. 1. The plot shows 0Jq — M c (solid line) and Md (dotted line), as functions of 
r in . The graphs are obtained for V = 5/3, and different values of M c = 10, 10 _1 
and 10 -3 respectively. 



Graphs of ujq — M c for different values of M c and different polytropic 
exponents T are shown on Figs. [3] and [4} We observe that this difference 
loq — M c , treated as a function of fj n , tends to zero for M c — > oo (test-fluid 
limit), and to a fixed function cD^fm, M c = 0) for M c — > 0. The solution with 
M c = can be obtained numerically, and tI>g(fi n ,M c = 0) has a maximum 
of the order of 0.13 for T = 5/3 and 0.06 for T = 4/3. That, in turn, yields 
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Fig. 2. The same as in Fig. [IJ but for V = 4/3. 



upper estimates for the ratio Cjq/M c , namely 



Mr. 



< 1 + 



and 



M, 



0.13 
0.06 



for r 



3' 
4 



< 1 + ^ for r 

M r 3 



(10) 



(11) 



For the reasons of clarity, the values of uiq/M c — 1 are also plotted in Figs. [H] 
and pi The maximum value of the rescaled mass of the disk M& is attained 
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0.14 




Fig. 3. Dependence of Qq — M c on f- m for different values of M c and T = 5/3. 



for M c = 0. It is of the order of unity for T = 5/3 and even smaller for 
r = 4/3. Thus, it is clear from Eqs. (10) and ( |11[ ) that the central mass 
M c influences the angular velocity u much stronger than Af<j, or, to put it 
differently, M c enters the expression for uj with a larger weight than M^. 

Numerically obtained solutions correspond to thickest disks when M c = 
0, and become infinitely thin for M c — > co, that is in the test-fluid limit. For 
a fixed value of f- m , the relative height of the disk H/ (r- m — r out ) (we assume 
that the disk extends from z = — H to z = H) decreases monotonically 
with M c . Thus, in principle, geometrical parameters such as the ratio of 
r in/ r out and the relative height of the disk should allow one to identify a 
particular solution and to obtain values of Coq/M c and M^/M c , provided 
that the polytropic exponent is known. 

Figures [7-10 show the dependence of uiq/M c on the mass ratio M&jM c . 
We focus on the physically interesting range of parameters < M&jM c < 1 
and < fi n < 1/2. Results obtained for T = 5/3 are shown in Figs. [7] 
and [8] the data for F = 4/3 are depicted in Figs. [9] and [Toj These graphs 
confirm quantitatively what was anticipated in Sec. [T] For small values of 
f; n (elongated disks) the motion is strictly Keplerian with luq ~ M c , even if 
the mass of the disk is comparable to the central mass. Faster rotation can 
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Fig. 4. The same as in Fig. [3| but for T = 4/3. 



occur for ~ M c only if in addition f- m is large, but even for f; n = 1/2 we 
5/3 and uJq « 1.4 for T = 4/3. In general these results 



get 



1.3 for r 



weakly depend on the polytropic index T. 



6. Summary 

We show explicitly, by analysing thousands of solutions of polytropic 
hydrodynamic models, that the self-gravity of a gaseous heavy disk orbiting 
a central mass can be compatible with the Keplerian motion. Nevertheless 
the masses of the two components can be resolved, if in addition to the 
rotation velocity two geometric pieces of information are given, namely the 
ratio of the innermost radius r m to the outermost disk radius r out and its 
relative maximal height H/(r out — r[ n ). 

How legitimate is the (often made) guess that rotation is strictly Ke- 
plerian, ~ GM C , which allows one to infer the central mass from the 
observed rotation law? A result of Hure et al. [11] (obtained for dust) sug- 
gests that the error implied by this identification should be proportional to 
the ratio of Md/M c , where is the mass of the disk. We show that this is 
not true for polytropes. Our numerical results suggest that polytropic Ke- 
plerian disks rotate with the angular velocity exceeding (GM c /r 3 ) 1//2 , but 
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Fig. 5. Dependence of Wq/M c — 1 on f- m for different values of M c and T = 5/3. 



smaller than (G(M C + M ( j)/r 3 ) 1 / 2 . The lower bound in these inequalities 
can be obtained analytically for special configurations, reported in Sec. 3. 

Presented results are based on hydrostationary modelling. Situation can 
be different for systems strongly influenced by physical factors that were 
not discussed here. For instance, models of strongly radiating thick toroids 
considered in [TO] are not only Keplerian, but strictly Keplerian with ojq = 
GM C , although those disks are massive, and the value of n n /r ou t is not small. 
In other words, very strong radiation can invalidate conclusions concerning 
hydrostationary systems. On the other hand, rotation velocity of weakly 
radiating disks analysed in |14j is consistent with the results of this paper. 

This analysis is entirely Newtonian, and therefore it should be repeated 
in the general-relativistic context. We do not expect essential changes in 
those systems, where the inner disk boundary is well separated from a black 
hole horizon, r; n 3> 6M C . This issue - and further related questions - de- 
serves separate studies. 
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Appendix A 

We derive in this appendix various criteria - restrictions on equations of 
state - that guarantee finite disk configurations. The configuration is called 
finite if it can be contained in a compact subset of M 3 . We redo the proof, 
since the existing versions [TB] do not include the case with the singular 
central potential —GM c /\x\. 

We start with an observation that the inifiniteness in the axial direction 
z is incompatible with the cylindrical Keplerian rotation law (any cylindrical 
rotation in fact). This can be demonstrated as follows. Assume that the fluid 
extends to infinity along z direction for two different constant radii r± and 
T2- Then from Eq. ^ we have 

lim (/t + $ + * c ) = $c(n) = C, 

z— >oo, r=n 

and the same applies to r = r%. Thus, we would conclude that ^(n) = 
^c{ r 2) = C, which contradicts the assumption that the motion is Keplerian. 
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Fig. 7. Dependence of u)q/M c on M d /M c for different values of f- m and T = 5/3. 

In principle the disk could still extend to infinity at the equatorial plane. 
Assume now that the disk has a finite mass. Then the only value of constant 
C appearing in Eq. Q allowing for an infinite configuration is C = (note 
that the Keplerian centrifugal potential is normalised to zero at r — > oo). 
This situation can be excluded for some cases with the help of the virial 
theorem ([7]). 

The first necessary step is to establish suitable function spaces for the 
density, and the pressure, that include infinitely extended fluids, and for 
which the virial theorem can still be proven. This can be done in terms 
of weighted Sobolev spaces. Details (which we omit) can be found in |15j . 
Apart from the suitable falloff behaviour we also assume that the support 
of the density is a connected set. 

Multiplying the integrated Euler equation by p and integrating over M 3 
we get 



C J pd 3 x = J phd 3 x + J p$ g d 3 x + / p 
Combining this equation with the virial theorem ([7]) one obtains 
C J pd 3 x = i J p$ g d 3 x + j(ph- 3p)d 3 x. 
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Fig. 8. Dependence of Wq/M c on Md/M c for different values of f- m and T = 5/3. 

Now, for those barotropic equations of state for which ph — 3p < every- 
where, we can conclude that C is strictly negative, and hence the fluid must 
be finite. For polytropic equations of state we have ph — 3p= ((3 — 2r) /(T — 
l))p, and the finiteness is guaranteed for T > 3/2. 

The condition on the equation of state can be relaxed on the expense 
of additional assumptions on the rotation. Making use of the virial theorem 
we can also write 



C J pd 3 x = j{ph~ 6p)d 3 x - J 



P 



GM C 



x 



d 6 x. 



If (jJq/t — GM c /\~x\ > everywhere, it is enough to demand that ph — 6p<0 
in order to conclude that the fluid must be finite. For polytropes the later 
condition simply states that T > 6/5. 



Appendix B 

Below we report results of the convergence tests for typical disk config- 
urations obtained with our SCF numerical scheme. In our implementation, 
numerical precision is controled by the grid resolution, the maximum num- 
ber of Legendre polynomials used in the angular expansion L, and a value 
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Fig. 9. The same as in Fig. [7j but for V = 4/3. 

of the maximal difference between density distributions obtained in the last 
two consecutive iterations ptoi- (In each iteration we compute the quantity 
/Jerr = niaxjj \pf^ —pfj\ - Here index k numbers subsequent iterations; in- 
dices i and j refer to different grid nodes. The iteration procedure is stopped, 
when perr < p tol .) 

Table [T] shows the dependence of results (ci)o and the virial parameter 
e v ) on the grid resolution, the maximum number of Legendre polynomials 
L, and p t oi- For this test we have assumed M c = 1, f- m = 0.3, T = 5/3. The 
corresponding mass of the disk yields M<± « 0.29. 

In addition, Tables [2] and [3] illustrate the dependence of the convergence 
properties on the value of M c (or, equivalently, the ratio of M c /M&). Data 
listed in these tables were obtained for L = 200, T = 5/3, r m = 0.3, and 
p tol = 10~ 5 . Table § corresponds to M c = (in this case M d w 0.49). In 
Table|3| we assumed M c = 10 (M d ps 1.4 x 10~ 2 ). 

Notice that the obtained value of the virial test parameter depends 
mainly on the resolution of the numerical grid. Values obtained in this paper 
agree with those reported by Axenov & Blinnikov [3 J , who were also testing 
their implementation on high-resolution grids. 
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Fig. 10. The same as in Fig. [sj but for T = 4/3. 
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Table 1. Typical dependence of the results on the grid resolution, the maximum 
number of the Legcndrc polynomials L, and the tolerance coefficient p to \- 



Resolution 


L 


Ptol 


W 


e v 




100 x 100 


200 


10" 


-b 


1.029945714 


2.22 x 10" 


-b 


200 x 200 


200 


10' 


-5 


1.029948775 


5.49 x 10" 


-6 


400 x 400 


200 


10" 


-5 


1.029949540 


1.30 x 10" 


-6 


800 x 800 


200 


10" 


-5 


1.029949076 


2.45 x 10" 


-7 


1600 x 1600 


200 


10" 


-5 


1.029949123 


1.67 x 10" 


-8 


3000 x 3000 


200 


10" 


-5 


1.029949140 


7.92 x 10" 


-8 


5000 x 5000 


200 


io- 


-5 


1.029949137 


9.51 x 10- 


-8 


100 x 100 


200 


10" 


-6 


1.029945714 


2.22 x 10" 


-5 


200 x 200 


200 


10" 


-6 


1.029948775 


5.49 x 10- 


-6 


400 x 400 


200 


10" 


-6 


1.029949540 


1.30 x 10- 


-6 


800 x 800 


200 


10" 


-6 


1.029949076 


2.45 x 10- 


-7 


1600 x 1600 


200 


10" 


-6 


1.029949123 


1.67 x 10- 


-8 


3000 x 3000 


200 


10" 


-6 


1.029948709 


2.54 x 10- 


-8 


5000 x 5000 


200 


10" 


-6 


1.029948705 


4.13 x 10" 


-8 


2400 x 2400 


100 


10" 


-b 


1.029948963 


7.35 x 10- 


-7 


2400 x 2400 


200 


io- 


-5 


1.029949132 


6.53 x 10- 


-7 


2400 x 2400 


300 


10" 


-5 


1.029949145 


6.47 x 10- 


-7 


2400 x 2400 


400 


10" 


-5 


1.029949147 


6.46 x 10- 


-7 


2400 x 2400 


500 


10" 


-5 


1.029949148 


6.45 x 10- 


-7 



Table 2. Dependence of the results on the resolution of the grid for a massive disks 



with M c 


= 0. Here p tol = 10" 5 


and L = 200. 


Resolution 


tl> 


e v 




100 x 


100 


0.339018059 


1.32 x 10- 


-4 


200 x 


200 


0.339031477 


3.67 x 10- 


-5 


400 x 


400 


0.339034828 


1.29 x 10- 


-5 


800 x 


800 


0.339035663 


6.99 x 10- 


-6 


1600 x 


1600 


0.339035431 


5.51 x 10" 


-6 


3000 x 


3000 


0.339035495 


5.15 x 10" 


-6 


5000 x 


5000 


0.339035486 


5.06 x 10- 


-6 
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Table 3. Dependence of the results on the resolution of the grid for a light disk 
with M c = 10. Here gtoi = 10~ 5 and L = 200. 
Resolution cDq e v 



100 


X 


100 


3.164121422 


1.14 x 10- 


-6 


200 


X 


200 


3.164121530 


2.32 x 10" 


-7 


400 


X 


400 


3.164121439 


1.25 x 10" 


-9 


800 x 800 


3.164121436 


5.66 x 10" 


-8 


1600 


X 


1600 


3.164121424 


7.12 x 10" 


-8 


3000 


X 


3000 


3.164121425 


7.46 x 10- 


-8 


5000 


X 


5000 


3.164121425 


7.55 x 10- 


-8 



